Zero-bias conductance in carbon nanotube quantum dots 
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We present numerical renormalization group calculations for the zero-bias conductance of quan- 
tum dots made from semiconducting carbon nanotubes. These explain and reproduce the thermal 
evolution of the conductance for different groups of orbitals, as the dot-lead tunnel coupling is var- 
ied and the system evolves from correlated Kondo behavior to more weakly correlated regimes. For 
integer fillings = 1, 2, 3 of an 5'!7(4) model, we find universal scaling behavior of the conductance 
that is distinct from the standard SU{2) universal conductance, and concurs quantitatively with 
experiment. Our results also agree qualitatively with experimental differential conductance maps. 
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Introduction: Carbon nanotubes have generated im- 
mense interest due to their rich transport properties [l|. 
Their small capacitance generates a Coulomb blockade 
regime at low temperatures jll, and a single electron 
transistor can be made out of weakly coupled nan- 
otubes. The Kondo effect, studied in transition metal 
ions and rare earth compounds for some fifty years 0, @] , 
is now equally a classic hallmark of many-body physics in 
nanoscale devices, including carbon nanotubes [61] as well 
as semiconducting quantum dots, molecules and mag- 
netic adatoms on metallic surfaces (for a review see j?! ) • 

While the 'standard' Kondo effect arising in nano- 
devices is the orbitally non-degenerate, spin-i {SU{2)) 
Kondo effect it is known that in quantum dots 
made from high quality nanotubes, a series of two spin- 
degenerate orbitals originates from the two electronic 
sub-bands [1, Consecutive filling of these orbitals 
should thus yield an SUiAUype Kondo effect d [13, El 
II, III Q 111, [S E H E . This has indeed been 
observed in recent experiments 
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including in 

particular systematic exploration of the SU (4) conduc- 
tance regimes via careful control of sample contact 
transparency. 

In this paper, we present a consistent picture which 
explains the experimental zero bias conductance [l2| . 
We study the thermal evolution of the conductance as a 
function of applied gate voltage, for an SU (4) Anderson 
model, using Wilson's numerical renormalization grour 



moaei, using Wilson s numerical renormalization group 
(NRG) ]^UM (previous NRG work [3, El, [3, [13, H, El 
has been confined to T = 0); and show that on progres- 
sive increase of the dot-lead tunnel couplings, the sys- 
tem exhibits a rich evolution from correlated Kondo be- 
havior with strong associated Coulomb blockade peaks 
in the zero-bias conductance, through to a weaker cou- 
pling regime where charge fluctuations are significant and 
Coulomb blockade oscillations suppressed. We obtain 
universal scaling functions for the zero-bias conductance, 
which deviate significantly from standard SU{2) Kondo 
scaling (reflecting the different universality class of our 
model); and show that this captures experiment quanti- 
tatively. 



Theory: Interacting quantum dots, molecular junc- 
tions and other nano-devices are modelled by the in- 
teracting region Tiimp, a set of non-interacting reser- 
voirs I-Lb and a coupling between the subsystems T-Lt'- 
H = TLimp + Ti-B + 'Ht- For the carbon nanotube 
quantum dot we restrict ourselves to the filling of one 
shell comprising two degenerate orbitals, and consider 
the local Hamiltonian TLimp — \U{N — Ng)'^. Here 

^ = J2aa'^iadca = Ea^^"'^' ^hcrc 4„ crcatcs a CT- 
spin electron in orbital a — 1,2, the charging energy 
U = e^/C = Ec with dot capacitance C, and Ng denotes 
the dimensionless external gate voltage. Equivalently, 



(1) 



with level energy Ed — ^{l — 2Ng) and m = (a, tr) a fla- 
vor index. For the isolated dot, with integer ground state 
charge {{N) =) Nd = — 4, the edges in the Coulomb 
blockade staircase between charge Nd and Nd+l occur for 
half-integral iVg = i(2Ar<j+l) {ie Ed = -UNd). InEq.(l) 
we have taken U' = U, with U [W] the intra- [inter-] 
orbital Coulomb repulsion, and have also neglected an 
exchange splitting, J. This is consistent with the argu- 
ment [24I that for carbon nanotubes, U' ~ U and J is 
much smaller than the charge fluctuation scale; and is 
known to be compatible with the nanotube data (23j . 

The tunneling to the two leads v = L,R\s assumed to 
be spin and orbital conserving 12, 3], 
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such that the overall H has S'[/(4) symmetry. Only the 
binding combination of lead states, Ckaa = 7 ^u^kaav 
with t = [^i +i|j]^, then couples to the dot. We can thus 
drop the lead index v and consider one effective lead: 



(3) 



Note however that the different couplings still enter via 
the conductance prefactor Go ~ {e^ /h)AT iX n/ l + 
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Tij)^, which reaches the unitary hmit of /h for a perfect 



symmetric channel = Ffl [2j|, where = irt^p with 
p the density of states of the leads at the Fermi energy. 

Method: We solve the Hamiltonian accurately using 
the powerful NRG approach j5|]; for a recent review see 
2l|. Finite temperature NRG Green's functions are cal- 
culated using the recently developed algorithm 25|, em- 
ploying a complete basis set of the Wilson chain. Origi- 
nally derived for real-time dynamics of quantum impuri- 
ties out of equilibrium , this ensures that the spectral 
sum-rule is exactly fulfilled and the NRG occupancy ac- 
curately reproduced. The self-energy of the Green func- 
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tion is calculated in the usual way 

The local Green function determining transport is 
given by G^(2:) = [z-e- A(z)-E(z)]~i with z = w+iO", 
A = Ylin it Sfe[-2^ ~ the dot-lead hybridization, and 

E(z) the interaction self-energy; the local spectrum is 
D(lu,T) = ^^mGi(z). At T = 0, the Friedel sum 
rule (271, |28| relates the number of impurity electrons 
Nimp to the conduction electron scattering phase shift 
S — irNimp/N, where N is the number of channels (here 



N 



4). Nimp = Nd + AN includes all electrons on the 
quantum dot {Nd = {J2a a ''^aa)) and the number of dis- 
placed electrons in the bath AN [2^ , and can be obtained 
directly from the NRG [5^. For a vanishing 3mEcrQ(«0~) 
as T — > 0, the zero bias conductance must approach 



lim G(T) = 4Gosm^(S) = 4Gosin2 



N 



■imp 



(4) 



which is indeed recovered in our calculations. 

Results: We use a symmetric, constant lead density 
of states po = ^/{'^'D) with a band width D = SOFq, Fq = 
ImeV ~ 11. being the energy unit used throughout. 
We take U — IOFq = IQmeV, corresponding to the ex- 
perimental estimate of the charging energy Ec = U [l3| ■ 

One of the striking, poorly understood features of the 
recent experiments on semiconducting carbon nanotube 
dots [l3| is the thermal evolution of the zero bias conduc- 
tance as a function of gate voltage. In the experiment, 
the orbitals are filled in groups of two (there being four 
such groups, I-IV [l3|)- In addition, the tunneling matrix 
elements increase with gate voltage. We restrict ourselves 
to a single representative group, keeping the tunneling 
matrix elements constant for simplicity; and calculate 
the zero-bias conductance as a function of the dimen- 
sionless gate voltage Ng, for different tunneling strengths 
F = Fl -I- Tfj and temperatures, using the NRG spectral 
functions: G{T)/Go = -4 J^^dco {df{uj)/duj)t{uj) with 
t(u;) = ttTD(uj) the t-matrix and /(w) the Fermi func- 
tion [23|. The results are given in Figs. [1] and [2] for three 
different couphng strengths F = O.SFq, Fq, 2Fo. 

Our results concur well with experiment. Those for 
F = 0.5meV (Fig. [T^) track the T-evolution of Group 



I depicted in Fig. 2 of Ref. [I4I. Clear, pronounced 
Coulomb blockade peaks and valleys are seen in Fig. [T^. 
These disappear only at temperatures much lower than 
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FIG. 1: (color online) Zero bias conductance vs dimensionless 
gate voltage Ng, for different T and two tunneling strengths 
F = 0.5mel/ ((a)), ImeV ((b)); with V = Ec ^ lOmeV. The 
dotted line gives the T = values, Eg. jH); Nimp is also shown 
(dashed line). NRG parameters A = 3, Ns= 2000. 



the experimental base T — 1.3K , indicative of strong cor- 
relations {U/r — 20 here) and hence small Kondo scales 
Tk in the middle of the valleys {Ng = 1 (3) and 2). The 
distance between the peaks is seen to become smaller 
with decreasing T, also in agreement with experiment. 
The origin is a pinning of the many-body Kondo reso- 
nance in D{uj) close to the chemical potential at low-T, 
evident in the t-matrix depicted in Fig. [5] (discussed be- 
low). For T = 0, the conductance from the Friedel sum 
rule Eq. ^ is shown as a dotted line in Fig. [1^. Abso- 
lute values of G{T) also agree well with experiment. For 
Group I 1^ the coupling asymmetry is Fl/F^ « 2.5, 
and hence Gq — 0.8e^//i. With this, eg the calculated 
A^g = 1.5 (or 2.5) Coulomb blockade peak in G(T) for 
T — 1.65Ar translates to Gpeak ~ 1.7e^//i, in good agree- 
ment with experiment; and the widths of the peaks also 
concur with experiment. The slight experimental asym- 
metry of the Ng — 0.5 and 3.5 Group I conductance 
peaks can be attributed to a small increase in the tunnel 
coupling within the Group on increasing the gate voltage. 

Resuhs for F = ImeV (C//F = 10, Fig^ [It) resem- 
ble the experimental Group III orbitals [ij]- At higher 
T, remnants of the Coulomb blockade peaks remain vis- 
ible; they are absent at base temperature. While for 
stronger correlations {eg Fig. [1^) the charge steps are 
quite pronounced, and the different orbital filling regimes 
well separated in the T = conductance, for U /V — 10 
orbital filling is almost continuous. The physical origin 
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FIG. 2; (color online) Zero bias conductance vs Ng for differ- 
ent temperatures and F — 2meV . The dotted line gives the 
T = values from Eq. (|4]). Other parameters as in Fig. [1] 



of such behavior is in part the 'blocking effect' known 
from the theory of multi-channel models in the context 
of rare earth materials in which Coulomb interac- 
tions lead to single-particle lifetime broadening of order 
NT, since N relaxation channels are present. In addi- 
tion, the SU{4) Kondo scale in the Kondo regime is much 
larger than for SU{2) as the number of channels enters 
the exponent, T^^^^^/D ~ [T^^^^^/Df/^ In con- 

sequence, an SU{4:) model enters the weakly correlated 
regime for larger U/T than an SU (2) symmetric model. 

Fig. [2] shows results for a more weakly correlated 
U/T = 5. This tracks the T-evolution of the Group IV 
orbitals Here, consistent with experiment, charge 

fluctuations are significant and Coulomb blockade peaks 
in consequence absent, with the G{T 0) limit Eqn. @ 
being reached in practice at relatively high temperatures. 

The T-dependence of G{T) for the middle of the 
Coulomb valleys, Ng = 2 and Ng = 1 (or 3), is shown in 
Fig. [3^,b for 5 values of F. Comparison of the two shows 
the characteristic energy scales in all valleys are of the 
same order (see also Fig. |4] inset): the Coulomb blockade 
valleys are thus filled simultaneously with decreasing T, 
as seen in experiment [l2j . This in turn provides strong 
evidence for U' = U, since only a very slight decrease of 
U' from the SU{4) point U' = U causes a sharp drop in 
the Ng = 2 Kondo scale towards SU{2) form lldl. while 
the Ng = 1 scale is by contrast barely affected [la. [l9j . 

Defining a characteristic low-energy scale Tq (oc Tk) by 



G{Tq)/G{0) = 0.5, Fig.[3t shows the Ng = 2 results from 
Fig. [li rescaled as G(T)/G(0) vs x = T/Tq. On progres- 
sively decreasing F the data collapse, over an ever in- 
creasing T/To-interval, onto a common curve g{x). This 
is the universal S'J7(4) scaling conductance for Ng = 2. 
That limit is well reached, for all T/Tq < 10^ shown in 
Fig. [3t, by F = 0.2meV; but even for the largest F we 
consider {2meV, U/T = 5), the conductance scales onto 
the universal g(x) over a significant x-range. Similar con- 
siderations arise for the Ng = 1 results, and the universal 
SU{4:) g{x) for Ng = 2 and 1 (or 3) are shown in Fig. H) 
For X > 1, g{x) for Ng — 2 and 1 asymptotically coincide, 
while for lower x the two are distinct (refiecting different 
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FIG. 3: (color online) T-dependence of G{T) for five different 
values of F and two different dimensionless gate voltages: (a) 
Ng ^2 and (b) iV,, = 1 (or iVg = 3). (c) Scaling the Ng = 2 
results, G(r)/G(0) vs x = T/Tq. 
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FIG. 4: (color online) S'f/(4) scaling conductance g{x) = 
G(r)/G(0) vs X ^ T/To - solid line for Ng = 2, dashed 
for Ng = 1 (or 3). The SU{2) scaling conductance is shown 
as a dotted line. Experimental results [T^] for Groups II, III 
and IV (with Ng = 2) are also shown. Inset: evolution of the 
low-energy scale To, for Ng — 2 (solid) and Ng = 1 (dashed). 



leading corrections to the 5*17(4) fixed point). The cor- 
responding universal g{x) for a symmetric SU{2) model 
is also shown. It is seen to have a quite different form, 
and is known jl2| not to fit well the carbon nanotube 
data. For the SU{A) model by contrast, our results agree 
well with experiment - as seen in Fig. 3] where data for 
Groups IL III and IV in the middle of the two-electron 
valleys [131 ^^e compared to the Ng = 2 SU{4) g{x). At 
sufficiently high T, departure from universal scaling must 
inevitably occur (as above, and just visible for the Group 
II, III data), but the data collapse remarkably well onto 
the scaling conductance [1^ . 

Finally, to assess approximately the effects of a finite 
source-drain bias Vsd = V , we neglect explicit depen- 
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FIG. 5: (color online) Color-coded contour plot in the {V, Ng)- 
plane of (a) ts{V,Ng) for T = 0.5meV at T = 2.8if, as in 
Fig. [Iji; (b) t(a; = eV) for F = ImeV as in Fig. [1]d. 

dence of the self-energy on V. Granted this, if one of the 
contacts is sufficiently open that it acts as a tunneling 
probe, the T = conductance G{V,Ng) oc t{uj = eV) 
with t{LL}) the t-matrix as above. This situation is ap- 
propriate to the Group III and IV orbitals of 12]. If by 
contrast the contacts are more symmetrically coupled, 
then G{VJ^g) cx ts{V,Ng) = ^[t{Lu = +yV) + t{uj = 
— iey)] [2J]. This we believe is the relevant case for the 
Group I orbitals of^. For F = 0.5meV (as in Fig.[T}i), 
Fig. [5^1 shows a contour plot of Gots{V,Ng) (calculated 
for T = 2.8K The striking similarity to the ex- 

perimental Group I differential conductance map (Fig. 2b 
of [l2|) is evident, including the clear Coulomb blockade 
diamonds. For F = ImeV (as in Fig.[l]D), Fig.[5lD shows a 
corresponding contour plot of G^tiuj = eV)] this in turn 
captures well the Group III conductance map [l2| . 

Conclusion: Motivated by recent experiments on car- 
bon nanotube quantum dots [l^ . we have studied the 
evolution of an SU{A) Anderson model from correlated 
Kondo to mixed valent behavior, on progressive increase 
of the dot-lead tunnel couplings. NRG results obtained 
are in compelling agreement with experiment J,2], from 
the thermal evolution and scaling behavior of the zero- 
bias conductance, to differential conductance maps; and 
show that an SU{A) Anderson model provides a remark- 
ably faithful description of carbon nanotube dots. 
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